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ABSTRACT 

We present results for a large number of gamma-ray burst (GRB) afterglow 
light curve calculations, done by combining high resolution two-dimensional rela- 
tivistic hydrodynamics simulations using RAM with a synchrotron radiation code. 
Results were obtained for jet energies, circumburst medium densities and jet an- 
gles typical for short and underluminous GRBs, different observer angles and 
observer frequencies from low radio (75 MHz) to X-ray (1.5 keV). We summarize 
the light curves through smooth power law fits with up to three breaks, covering 
jet breaks for small observer angles, the rising phase for large observer angles 
and the rise and decay of the counterjet. All light curve data are publicly avail- 
able via [http : // cosmo . nyu . edu/ af terglowlibrary[ The data can be used for 
model fits to observational data and as an aid for predicting observations by 
future telescopes such as LOFAR or SKA and will benefit the study of neutron 
star mergers using different channels, such as gravitational wave observations 
with LIGO or Virgo. 

For small observer angles, we find jet break times that vary significantly 
between frequencies, with the break time in the radio substantially postponed. 
Increasing the observer angle also postpones the measured jet break time. The 
rising phase of the light curve for large observer angle has a complex shape that 
can not always be summarized by a simple power law. Except for very large 
observer angles, the counter jet is a distinct feature in the hght curve, although 
in practice the signal will be exceedingly difficult to observe by then. 



1. Introduction 



Short gamma-ray bursts (SGRBs) a re likely produced b y neutron star-neutron star or 
neutron star-black hole mergers (see e.g. lEichler et al.lll989l for an early exploration of this 
idea). This makes them physically different from long duration GRBs, which result from the 
stellar collapse of a massive star. The distribution of GRB durations is therefore expected 
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(and found) to be bimodal rather than continuous (IKouvehotou et aLlll993[ ). Nevertheless, 
theoretical models of both types of GRB share many similarities, the most important of 
them being the formation of an ultrarelativistic jet. SGRB jet models generally probe a 
different part of the possible parameter space for such jets than long GRBs, although some 



overl ap exists with underluminous instances of the latter (e.g. GRB 100316D. IStarling et al. 



2OIII ). The overall energy release for SGRBs is of the order io^^~^° ergs (rather than 10^^ 
ergs), the circumburst particle densities of order 10~^ — 1 cm~^ (rather than 1 cm~^) and 
they are less collimated (although there is currently li ttle observational confirmatio n of the 
latter, in hydrodynamical models for long GRBs, e.g. iMacFadyen fc Woosleylll999l . the jet 
becomes collimated by passing through a dense stellar interior . This mechanism is absent 



for SGRBs). Reviews of SGRB science ca n be found inlNakad (120071) . iGehrels et al.l ((20091) 
and a recent comparison to long GRBs in iNysewander et al.l ( 120091 ). 



Like long GRBs, short and underluminous GRBs produce afterglows peaking at pro- 
gressively longer wavelengths with time, although they are harder to detect because they are 
intrinsically fainter. Analytical models of SGRB afterglows suffer from the same simplifica- 
tions and shortcomings as those of long GRBs (mainly with respect to jet decollimation and 
off-axis emission). These can be addressed through combining high-resolution relativistic hy- 
drodynamics (RHD) simulations with numerical radiative transfer for s ynchrotron radiation. 



Such simulations have already been performed for long GRBs (e.g. iZhang &: MacFadyen 
2OO9I : IVan Eerten et al.l l2010al ). A more accurate understanding of short GRB afterglows 
is currently especially interesting not only because they are actually being detected start- 
ing a few years ago, but also because a new generation of extremely sensitive detectors of 
SGRBs is becoming operational. On the one hand there are instruments such as LOFAR or 
SKA, that will detect SGRB (afterglow) emission through the traditional electromagnetic 
(EM) channel, but at unprecedented long wavelengths on the order of tens of MHz rather 
than GHz (The lower limit goa l for SKA is 60 MHz, for SKA pioneer project A SKAP it is 



300 MHz, I Johnston et al.ll2008l for LOFAR it is ~ 10 MHz, iRottgerind 120061 ). Transient 
monitoring campaigns with these instruments should be able to detect afterglows even if the 
prompt emission remains unseen. On the other hand, completely new channels are becoming 
available for GRB dete ction: multiple grav itational-w ave (GW) detectors a re currently in 
operation (e.g. LIGO, lAbbott et al.l l2009l and Virgo, lAcernese et al.l l2008l ) and upgrades 
are anticipated. The amount of information that can be obtained from GW detections can 
be significantly enhanced by in formation from their EM counterparts that can help break 
degeneracies in GW model fits (INissanke et al.ll2010l ). Also, the expected observer time be- 
tween the GW signal and the peak of the afterglow signal is expected to be on the order 
of several days at least, so a GW localization can be used to increase the odds of detecting 
an afterglow. It is therefore important to accurately understand the relationships between 
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SGRB energy, collim ation and observer angle and their observational implications (see also 



Nakar fc Piranll201l[ ) 



To further our understanding of SGRB afterglows we have performed a series of two- 
dimensional RHD simulations of SGRB jets interacting with the circumburst environment. 
In this letter we present, for the first time, short GRB afterglow light curves for observers 
positioned both on the jet axis and off-axis that are calculated from RHD simulations and 
include both synchrotron emission and synchrotron self-absorption. In section |2] we explain 
our methods. In section [3] we discuss light curves and spectra for a number of cases. All data 
are publicly available at http : / / cosmo . nyu . edu/ af terglowlibrary , but in addition to this 
we provide tables showing results for smoothly broken power law fits to the light curves for 
the cases under discussion. The fits should be taken as a brief summary of the overall 
shapes of the light curves, rather than a definitive description that completely captures the 
underlying physics. In section H] we draw conclusions. 



Simulations and radiation 



We have performed 12 RHD sim ulations in 2D spherical coo rdinates using the parallel 
RAM adaptive- mesh refinement code (jZhang fc MacFadyenI 120061 ) . four of which we will de- 
scribe in detail in this paper (the resulting light curves for the others are publicly available 
on the website). Taking a conic section from the Blandf ord-McKee (BM) self-simila r solution 
for a relativistic explosion in a homogeneous medium (IBlandford &: McKeelll976l ) as start- 
ing point, these simulations cover all possible combinations of the following values for the 
physical quantities determining the dynamics of the system: energy in both jets Ej = 10^^ 
or 10^° ergs, circumburst number density n = 10~^, 10~^, 1 cm~^, jet half opening angle 
6j = 0.2, 0.4 rad. Numerical resu lts for the dynamics o f sprea ding and decelerating BM jets 



have been described in detail in IZhang fc MacFadyenI (120091 ). Based on these findings, we 
have chosen a starting BM jet fluid Lorentz factor directly behind the shock front of 7 = 10 
(well in excess of l/0j, both for 6j = 0.2 and 0.4, when lateral spreading is expected to set 
in), a lab frame stopping time t = lOtNR and corresponding grid size r = ct (where c the 
speed of light and t^R as defined below) and a grid resolution such that the initial blast wave 
width R/T^ (with R the blast wave radius and F the shock Lorentz factor) was resolved by 
approximately 100 grid cells. In practice, the latter requirement lead us to 24 base level 
blocks (of 16 cells each) in the radial direction and 12 initial levels of refinement. We have 
used 2 base level blocks in the angular direction. As in earlier work, we gradually decreased 
the peak refinement level over time. In addition, we kept the peak refinement level for the 
inner regions of the jet a few levels lower than that of the outer regions. This avoids spending 
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too much computational effort on resolving Kelvin-Helmholtz type instabilities in the ffow 
that have little effect on the radiation and overa ll dynamics. The time t^R reffects the time 
when the blast wave was analytically expected Jpiran 2005 ) t o become n onrelativistic and 



settle into the self-similar Sedov- Taylor solution (jTaylorlll950t ISedovlll959l ) and is given by 



t 



NR 



-1/3 



days. 



Here Eiso,53 denotes isotropic equivalent energy in units of 10 ergs, related to the energy 
in both jets according to Eiso = 2Ej/9'j. The theoretical value for t^R was numerically 
fou nd to underestimate t h e tra nsition time and duration to spherical nonrelativistic flow 
and IZhang fc MacFadyeru (120091 ) found that the transition time is better approximated by 



^ 5tNR- 

We have combined simulation output (3000 snapshots per simulation) with a linear ra- 
diative transfer code that calculates the observed flux at various observer angles from rays 
through the evolving fluid, including synchrotron emission and synchrotron self- abso r pt ion . 
The current approach gen eralizes the radiation code described in I Van Eerten fc Wijerd (120091 ) : 
Van Eerten et al.l (j2010al ) to off-axis observers. It further differs from these studies in that it 



follow s the simplified approach to the general synchrotron emission used in Ivan Eerten et al. 



(120101 ) that uses a globa l approach to ele ctron cooling rather than a local one and that in 
turn has been based on ISari et al.l (119981 ). The total number of rays required for a single 
observation is calculated through a procedure analogous to adaptive mesh refinement, where 
each possible refinement doubles the number of rays in a single group (or "block" , containing 
16 rays) either in the r or direction, where r and (p are polar coordinates on the plane 
perpendicular to the direction of the rays (i.e. towards the observer). A local total of 9 
refinements (in any combination of radial and angular refinements), starting with 24 base 
level blocks in the radial direction (r = - fluid grid maximum) and 2 in the angular di- 
rection (0 = — tt) was found to be sufficient to converge on a fixed flux value. Although 
our method provides us with spatially resolved images for off-axis observers as well as fluxes, 
these will be presented elsewhere. 

As usual in afterglow modeling, a number of parameters is used to capture the radiation 
physics. The accelerated electron power law slope has been set to p = 2.5, the fraction of 
accelerated electrons C,n = 1-0, the energy in these electrons as fraction of the thermal energy 
eE = 0.1, the fraction of thermal energy in the magnetic field = 0.1. 

In order to ensure complete coverage at early observer times we have used an analytical 
implementation of the BM solution at fluid Lorentz factors > 10 rather than simulation 
output to calculate emission and absorption. Early time contributions have been confirmed 
to connect smoothly to those from simulation output. However, even with additional early 
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time contribution, the earliest observer time with effectively full coverage still differs between 
observer angles. 



3. Light curves and spectra 

We haved calculated light curves at the following four frequencies: 75 MHz (radio, 
LOFAR, SKA), 1.43 GHz (radio, WSRT, VLA), 4.56 x 10^^ Hz (R-band, VLT) and 
3.63 X 10^'' Hz (1.5 KeV X-rays, Swift XRT). In this letter we summarize and discuss the 
results in detail for the following cases: 



• Ei 



• E, 



10^^ ergs, dj = 0.2 rad, n = lO'^ cm^^ (A) 

10^^ ergs, Oj = 0.4 rad, n = IQ-^ cm^^ (B) 

10^° ergs, Oj = 0.4 rad, n = 10"^ cm"^ (C) 
Iq50 gj,gg^ Q, ^ g_4 j^g^^^ n = 1 cm"^ (D) 



This way we cover both small and large opening angles and the effect of increased jet energy 
and circumburst density. For all cases we have computed light curves for a range of observer 
angles: Oobs = 0, 0.1, 0.2, 0.4, 0.8 and tt/2 rad. 

The resulting light curves for case B have been plotted in Fig. [T] For three observer 
times we have calculated broadband spectra as well, and these are plotted in Fig. [2l For the 
whole set of simulations (case A-D), we have summarized the shapes of the light curves by 
fitting smoothly connected power laws in time, using 



t 



-soi/3o 



,^01 



si2(/3i-/32)' 



^23 



X 



S23(/32-ft) 



1 
=23 



(2) 



for all observer angles except 7r/2, when both jets are seen exactly on edge and 



t 



-soi/3o 



+ 



t 



-soi/3i 



(3) 



is sufficient. Using multiple smoothly connected po wer laws to describe the data is common 
both in theoretical and observational studies (e.g. iBeuermann et al.l Il999l : iGranot fc Sari 
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20021 ). Fit parameter Fq sets the scale of the hght curve (here, we have set redshift z = 
and observer luminosity distance = 10^^ cm). Different power law regimes meet at 
transition points toi; ^12 and ^23, measured in days. The slopes at the different regimes are 
given by /3o, /9i and (32, while the sharpnesses for the transitions are given by sqi, S12, S23. 

Fit results for case A-D are listed in tables [1] and El Complete datasets are publicly 
available (at http : / / cosmo . nyu . edu/ af terglowlibraryj) and we emphasize that the power 
law fits are meant only as convenient summary, rather than as a full description based on 
the underlying physics of the afterglows. Nevertheless, the first temporal break tgi can be 
roughly interpreted as the jet break time for small observer angles. For high observer angles 
it marks the differe nce between the ris e (where relativistic beaming dominates the shape of 



the light curve, e.g. iGranot et al.ll2002[ ) and decay of the signal. The second and third break 



ti2, ^23 are used to summarize the rise and decline of the counterjet. 

Fitting for 11 fit parameters allows for a lot of freedom, and in most cases the resulting 
fit function captures the simulation light curve with only an occasional difference of up to a 
few percent. In Fig. Owe have plotted a number of fit results for case B, that illustrate the 
accuracy of the power law fits. The fit results have been obtained using a straighforward 
implementation of nonlinear least squares fitting while assuming a fractional error of 10 
percent on the simulation points. The number of datapoints varies between ~ 200 for on- 
axis observers and ~ 100 for observers completely off axis. The number of degrees of freedom 
(dof) varies similarly, with 11 fit parameters (or 5, for a single break) subtracted from the 
number of datapoints. All fits had x^/dof < 1, except for a number of separately marked 
cases. Sometimes fits were insensitive to some of the fit parameters. Where this was the 
case (when the error on the fit parameter was comparable in size to the parameter itself), 
only a single digit is given in the table. 

From Figs. [T], [2] and the tabulated results, we draw conclusions regarding the general 
structure of SGRB light curves. For small observer angles, the r esults clearly confirm that 



afterglow jet breaks are chromatic, which was first reported in IVan Eerten et al.l (l2010bf ) 
(albeit there for long GRB's, both for ID top hat jet simulations and a medium-resolution 
2D simulation using a different hydrodynamics code). At low (radio) frequencies, the jet 
break is consistently postponed with respect to the jet break in the optical and X-ray. 
Although the fitted break times also differ between optical and X-ray, the current approach 
that approximates the electron cooling time by the explosion duration is not sufficiently 
realistic to allow for definitive statements on the chromaticity of the break between the two 
frequencies (a detailed treatment of electron cooling lies beyond the scope of this work. The 
practical relevance of off-axis X-ray afterglow light curves for SGRBs is limited). Although 
we have not accounted for this in our fits, the jet break splits into two breaks when the 
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observer moves off-axis. The Oobs = 0.2 rad light curves from Fig. [T] provide an example of 
this: at early times they have again joined their repective on-axis light curves. For observers 
exactly on edge, only the late time break is left, with the early time break having moved to 
t = 0. This wide range of jet breaks is unfortunate from a practical point of view. It means 
that an equation for the break time tj such as 



tj = 3.5(1 



^1/3 -1/3 



+ Gobs 

0.2 



8/3 



days. 



(4) 



from Ivan Eerten et al.l ( l2010l ) will be consistent with the data (as can be seen from a com- 
parison to values in tables [Hand [2]), but only given the wide range of jet break times across 
the frequencies. Nevertheless, the scalings in eq. SJare confirmed. 

All lig ht curves confirni the clear rise of the counterjet that was both ana lytically ex- 



pecte d (e.g. Granot &: Loeb 2003 ) and a robust feature of earlier numerical work (jzhang fc MacFadyen 



20091 ) . They show that at late times and for off-axis observers, synchrotron self-absorption 
(s.s.a.) has little effect on the light curves. For case B, this can be seen directly in Fig. |2l 
For small opening angles and radio frequencies, the influences of s.s.a. and the jet opening 
angle on the light curve can become hard to disentangle. The most extreme case is provided 
by case D, where the circumburst density n is the highest. This is not unexpected given that 
the synchrotron break frequency scales as u^a oc n^^^ when t/.^n , lies b elow the synchrotron 
break frequency Um and as oc n"^^^^ otherwise (IGranot fc Sarill2002l ). 



Light curves for large observer angles exhibit a steep rising phase that will be more 
complex than a straight power law, especially at low frequencies. The low radio light curve 
plots in Fig [1] for 9obs equal t o 0.8 and 7r/2 provide e xamples. This was not reported in 



earlier numerical studies (e.g. Ivan Eerten et al.l l2010l : iGranot et al.l l200ll . where the early 



rising part has been truncated from the light curve). It is a genuine feature that depends 
on the s pectral shape rather than lateral spreading and is also seen in simplified analytical 
models ( 



Rhoadslll999l : ISari et al.lll998l : IGranot et al.l 120021 ) I 



A long term feature common to all light curves not fully captured by the power law fits 
is the gradual transition from relativistic to nonrelativistic flow. At the last observer time 
covered the simulations have not yet been nonrelativistic sufficiently long for the expected 
nonrelativistic slope above of P = —1.65 to become dominant (note that at any given 
observer time the observed flux is the combined signal from a range of emission times). In 
practice, observing a SGRB afterglow fully in the nonrelativistic regime will be exceedingly 
unlikely. 
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4. Discussion and conclusions 

We have performed a series of high-resolution RHD simulations in 2D to calculate the jet 
outflow for physical parameters typical of those expected for subenergetic GRB's. From these 
we have calculated afterglow light curves at various frequencies, covering low radio (75 MHz) 
up to X-ray (1.5 keV) and for observer angles from to 7r/2 rad. The data for all light curves 



from this paper are publicly available via http://cosmo.nyu.edu/afterglowlibrary, that 



also provides results from a more extensive probe of parameter space. We summarize the 
light curves via smooth power law fits that capture features such as the jet break for small 
observer angles, the early time rise due to relativistic beaming for high observer angles and 
the rise and decay of the counterjet. The results here present the most accurate calculations 
to date of light curve predictions of the standard afterglow jet theory as it applies to short 
GRB's, fully accounting for aspects such as jet spreading, observer position and arrival time 
effects. Although we do not discuss this in detail in this work, the light curves show that 
SGRB / underluminous afterglows should in principle be observable (at least in the radio) 
even for observers outside the jet cone. The light curves in this paper and in the on-line 
database should prove useful for detectability estimates using future radio telescopes such 
as SKA and LOFAR. Such estimates will also benefit the gravitational waves community, 
since the amount of information that can be extracted from GW measurements increases 
significantly when EM counterparts are observed as well. Furthermore, GW observations 
can aid the search for EM counterparts. 

In general, the afterglow light curves are well described by smooth power law fits with up 
to three breaks, although sometimes the rising phase for high observer angles is problematic. 
The effects of increasing circumburst density and jet energy are as expected from theoretical 
models. The jet break for small observ er angles varies greatly between frequencies, confirming 



a result from lVan Eerten et al.l (j2010bl ). Increasing the observer angle postpones the jet break 



(it really splits the break into two separate breaks, but the second break is the strongest). 
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Table 1: Power law fit results for case A (top) & B (bottom). The observer angle is in radians. Flux level 
_Fo is in mJy. The break times tm, ti2 and ^23 are in days. The occasional lower case letters in the first 
column mark the following: a) Poor fit (meaning x^/dof > 1, with value up to ^ 5), caused by the complex 
shape of the initial rise of the light curve, b) Poor fit because the synchrotron break frequency Vm passed 
through the observed band around ten days and the corresponding complication for the light curve was not 
taken into account (see also Figs. [Hand [2]). c) Good fit, albeit with the break times ti2 and ^23 swapped 
and the final slope rising, hindering an interpretation of (5^ as the final slope. When a fit was insensitive to a 
fit parameter and the resulting error on the fit parameter of the same order as the parameter itself, a single 
digit has been used. All light curves for case A run until 4,000 days observer time, with starting observer 
times of 1 day (0.0, 0.2, 0.4 rad), 15 days (0.4 rad) and 50 days (7r/2 rad). For case B, the latest observer 
time covered is 3,200 days, with starting times of 1 day (0.0, 0.2, 0.4 rad), 6 days (0.4 rad) and 40 days (7r/2 
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Table 2: Power law fit results for case C (top) & D (bottom). The observer angle is in radians. Flux level 
Fq is in ni.Jy. The break times ioi, ti2 and arc in days. Lower ease letters in first column mark the 
following: a) Poor fit (meaning x^/dof > 1), caused by the complex shape of the initial rise of the light 
curve, d) Good fit, but the first two breaks represent complexities due to the combined effects of the jet 
structure and self-absorption, rather than jet break and rise of the counterjet. e) Good fit, but the double 
peak feature from the combined effect of jet structure and self- absorption is missed and the break times toi 
has no straightforward interpretation, f) There was such a clear break here in the rising phase that we fitted 
the curve in two parts: 15-600 days and 300-1500 days. When a fit was insensitive to a fit parameter and 
the resulting error on the fit parameter of the same order as the parameter itself, a single digit has been 
used. All light curves for case C run until 10,000 days observer time, with starting observer times of 2 days 
(0.0, 0.2, 0.4 rad), 30 days (0.4 rad) and 150 days (7r/2 rad). For case D, the latest observer time covered is 
1,500 days, with starting times of 1 day (0.0, 0.2, 0.4 rad), 3 days (0.4 rad) and 15 days (7r/2 rad). 
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Fig. 1. — Observed luminosity light curves for Ej = 10^^ ergs, 6j = 0.4 rad, n = 10"^ cm~^ 
(case B). Observer frequencies from top to bottom: 3.63 x 10^^ Hz, 4.56 x 10^^ Hz, 1.43 GHz 
and 75 MHz. The legend applies to all plots. 10 days, 50 days and 1 yr have been marked 
with vertical dotted grey lines. Spectra for these times are provided in Fig. [2J 
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Fig. 2. — Spectra for Ej = 10^^ ergs, 6j = 0.2 rad, n = 10 ^ cm ^ at tots = 10 days, 50 days 
and 1 yr (top to bottom plot), for various observer angles. The legend applies to all plots. 
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Fig. 3. — Direct comparison between power law fits and simulated light curves (broad light 
grey curves) for case B. The legend applies to both plots. The top plot shows optical 
(4.56 • 10^^ Hz) light curves, the bottom plot shows low radio (75 ■ 10^ Hz) light curves. 



